Entropy production, energy dissipation and violation of Onsager relations 

in the steady glassy state. 
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In a glassy system different degrees of freedom have well-separated characteristic times, and are 
described by different temperatures. The stationary state is essentially non-equilibrium. A gener- 
alized statistical thermodynamics is constructed and a universal variational principle is proposed. 
Entropy production and energy dissipation occur at a constant rate; there exists a universal relation 
between them, valid to leading order in the small ratio of the characteristic times. Energy dissipation 
(unlike entropy production) is closely connected to the fluctuations of the slow degree. Corrections 
due to a finite ratio of the times are obtained. Onsager relations in the context of heat transfer are 
also considered. They are always broken in glassy states, except close to equilibrium. 
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Statistical thermodynamics is a universal and powerful 
theory for describing equilibrium states. It was general- 
ized to weakly non-equilibrium states, in an approach 
first started by Onsager, and further developed exten- 
sively, see e.g. jlj. 

It was recognized long time ago that concepts and 
methods of statistical thermodynamics can also be ap- 
plied to glassy non-equilibrium states Q. In such sys- 
tems the relaxation times depend strongly on tempera- 
ture. When cooling at a proper rate (varying from 10~ 2 
K/s for window glass to 10 5 K/s for metallic glasses) the 
equilibrium relaxation time becomes very large near the 
experimentally denned glassy temperature T g . The thus 
reached metastable state is not in equilibrium but, never- 
theless, can be described by a generalized thermodynam- 
ics, assigning different temperatures (so-called effective 
or Active temperatures) to processes with well-separated 
characteristic times Q - In spite of much progress 
in this area many important questions are still not fully 
understood. In particular, it concerns dissipative effects. 
However, as the steady glassy state is non-equilibrium, 
there exist a constant-rate entropy production, an energy 
dissipation, and a transfer of heat. Although the physi- 
cal importance of entropy production was stressed in the 
fundamental review Q , its systematic investigation has 
been continued only very recently |J. 

We shall consider the steady non-equilibrium glassy 
state of systems in which the subsystems are coupled 
to baths at different temperatures. The glassiness here 
is solely a consequence of the assumed large separation 
of time-scales of the subsystems. Our purposes are the 
following, (i) Derive the glassy stationary statistical dis- 
tribution and the corresponding thermodynamics, (ii) 
Propose a general variational principle for glassy ther- 
modynamics. (Hi) Investigate the entropy production 
and energy dissipation in the stationary non-equilibrium 
glassy state, (iv) Show the breakdown of the Onsager 



relations for heat transfer. 

To adstruct our conclusions, let us introduce the sim- 
plest glassy system. It consist of a pair of coupled 
stochastic variables x\, xi with Hamiltonian H{x\,x<i), 
which interact with different thermal baths and have dif- 
ferent characteristic time-scales. Such an approach pre- 
tends to establish all important, necessary ingredients 
of glassy behavior. A theory of statistical systems in- 
teracting with different thermal baths was investigated 
in jjj. The essentially new points of our approach are 
the large separation between characteristic times, and 
arbitrary difference between temperatures. The (over- 
damped) Langevin equations for the dynamics read: 



Ti±i = -diH 
1,3 = 1,2 
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where Ti, T 2 are the damping constants, and di — d/dxi- 
The Einstein relation between the strength of noise and 
the damping constant holds in Eq. ([!]) because the ther- 
mal baths themselves are in equilibrium |2j. It is as- 
sumed that the relaxation time toward the total equilib- 
rium (where T 2 = Ti) is much larger than all considered 
times; thus for our purposes T 2 and T\ are constants. 

Hereafter we shall assume that x 2 is changing much 
more slowly than x\\ this condition is ensured by the 
condition 7 = T1/T2 <C 1. Let us first indicate how 
the stationary distribution can be obtained to order 7 , 
which will give us the basic formulation of the general- 
ized glassy thermodynamics. Eqs. (|l|) can be investi- 
gated by the method of adiabatic elimination (Born- 
Oppenheimer method). First Eq. (|l]) for x\ is solved 
keeping the x 2 fixed, valid on relatively short time-scales 
where only Eq. (|I]) for x\ is relevant. In this case the 
Langevin equation has the equilibrium distribution 



P (xi\x 2 ) 



Z{x 2 ) 



exp(-/3 1 H(x 1 ,x 2 )), 



(2) 
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where Z(x 2 ) is the partition sum for a fixed value of x 2 . 
X\ can be carried out. At quasi-equilibrium of the x 2 - 
subsystem this average should be performed using the 
distribution (0). In this way we get from Eqs. (|l|) 
a related dynamics for the slow variable, in which the 
two particle Hamiltonian H(x±, X2) is replaced by the ef- 
fective one-particle Hamiltonian — Tiln Z(x 2 ). We thus 
have the effective equation of motion 



r 2 ±2 



_d_ 

dx' 



T 1 \nZ(x 2 ) + r l2 (t) 



(3) 



As the noise is due to a bath at temperature T 2 , see Eq. 
(flh, the equilibrium distribution of this equation reads 



Po(x 2 ) 



Z T ^ T -{x 2 ) 



Z = J dx 2 Z T ^ T ^(x 2 ) 



(4) 



The joint distribution of x\ and x 2 can now be writ- 
ten as Pq(xi, x 2 ) — Pq(xi\x 2 ) Po(x 2 ) . A similar approach 
is applied in spin-glasses and other disordered systems 
where n — T\/T 2 is considered as "dynamically gener- 
ated" replica number p"l| ]. Keeping this in mind we now 
consider the general situation. 

(i) If the state of a system is described by a distri- 
bution P(xi,x 2 ) — P{xi\x 2 )P{x 2 ) then there exist the 
general definition for the mean energy and entropy [Q: 
U = (H), S = -(InP). This latter Boltzmann-Gibbs- 
Shannon formula corresponds to the general statistical 
definition of entropy, relevant also outside of equilibrium 
[0]. The total entropy decomposes as S = Si + S 2 , where 

51 = J dx 2 P(x 2 )[- J dxiP(x 1 \x i )ha.P(xi\x 2 )], 

5 2 = - J dx 2 P(x 2 )lnP(x 2 ). (5) 

Si is the entropy of the fast variable X\ , averaged over 
the quenched slow variable x 2 , and S 2 is the entropy of 
the slow variable itself. This general result of the statis- 
tical thermodynamics can be again applied in our case 
when P(xi,x 2 ) = P (x u x 2 ). From Eqs. (|), (|), §) 
an important relation can be derived which generalizes 
the usual thermodynamical relation for the free energy 
F = -T 2 \nZ: 



F = U- T x 5i - T 2 S 2 



(0) 



This agrees with the expression of the free energy for a 
glassy system put forward previously by one of us . 
In that approach the equivalent of T 2 is the dynamically 
generated effective temperature, while here it is the tem- 
perature of a bath. The first and second laws of thermo- 
dynamics take the form 



dU =dQ+dW < TidS 1 ! + T 2 dS 2 +dW, 



(7) 



where dW is the work which is done on the system by 
external forces, and the equality in Eq. (^) is realized 



for a reversible process. Eqs. (|[^) are the manifesta- 
tion of the glassy thermodynamics which generalize the 
usual one to the case of non- equilibrium systems with 
well-separated time-scales. Here they have been obtained 
from Langevin equations under the sole assumption of a 
separation of time scales. 

(ii) Let us indicate how a variational principle can 
be obtained from a more general consideration. The 
usual Gibbs distribution for homogeneous equilibrium 
states can be obtained either from maximizing the en- 
tropy, keeping energy fixed, or from minimizing the en- 
ergy, keeping the entropy fixed. For the glassy state, 
which is non-homogeneous and out of equilibrium, one 
can minimize the mean energy, keeping both entropies 
Si and S 2 fixed (somewhat similar to the microcanoni- 
cal approach) . Following the standard method we should 
minimize, with respect to P{x 2 ) and P(xi\x 2 ), the La- 
grange function 

C = J dxidx 2 P(xi,x 2 )H + T 2 J dx 2 P(x 2 )\nP(x 2 ) 

+Ti ( dx 2 P{x 2 ) ( dxiP(a;i|a; 2 )lnF(xi|x 2 ), (8) 



where Ti and T 2 are Lagrange multipliers, and normalize 
the solutions. We then recover Eqs. (||), (|]) but now 
on the basis of more general variational principle. This 
clearly demonstrates the conceptual difference compared 
to the usual (local-equillibrium) thermodynamics. 

( in) Due to a difference between Ti and T 2 there is con- 
stant current of heat through the system. This implies a 
constant production of entropy and dissipation of energy. 
We investigate these effects taking into account possible 
7 corrections. The Fokkcr-Planck equation which cor- 
responds to Eqs. (0) reads 

2 

d t P(xi,x 2 ;t) + '^2d i J i = 0, 

i=l 

T i J i = P( Xl ,x 2 ;t)diH + T i d i P(x 1 ,x 2 ;t) (9) 

where Ji , J 2 are the currents of probability. The station- 
ary probability distribution can be expressed as 

Pi{x u x 2 ) = P {x u x 2 )(l ~ jAixux^) + 0{j 2 ), (10) 

The boundary conditions are, as usual, that P{xi,x 2 ) 
and its derivatives vanish at infinity. A is obtained from 
the stationarity condition dtP = 0, taking into account 
the orthogonality condition J dxidx 2 AP^ = 0, and con- 
sistency with 0(j 2 ) terms. The general expression for A 
is rather lengthy, but for a concrete model it is given in 
Eq. (|l9|). In the first order of 7 the steady currents are 
given by 



Ji^j^-PodiA, 
1 1 



■h=l%^Po5F 2 , (11) 
i ii 1 
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Notice that for J 2 the object A is not needed, but only 



SF 2 = -d- 2 H 



dyPo(y\x2)d 2 H(y, x 2 ), 



(12) 



being the difference between the force acting on the sec- 
ond subsystem and its conventional mean value obtained 
by averaging over the fast degree of freedom. Therefore 
some further results can be derived without knowledge 
of A, though it is needed for consistency checks and j 2 - 
corrections. The change of total entropy reads 

dStot = ds +as M +as 6 , 2 = dS- ftdiQ - /? 2 a 2 Q (13) 

where S = S\ + S2 is the entropy of the system defined 
by @, Sb,i, Sb t 2 are the entropies of the correspond- 
ing thermal baths, and diQ, d 2 Q are the amounts of 
heat obtained by the system from the thermal baths. 
Of course, from the conservation of energy we have 
tiiQ +&Qbath,i — 0, while dQb.i = Ti dSb.i holds because 
the baths are in equilibrium. Further, the expression 



diQ 

dt 



J dxidx 2 H(xi,X2)diJi 



(14) 



can be obtained from (^]). The entropy and the mean 
energy of the stationary state are constant: S = 0, 
Qi + Q2 = 0. Nevertheless, there exists a constant- 
rate transfer of entropy to the outside world (the thermal 
baths), and a stationary flux of heat through the system: 



Stat = (fil ~ 



(15) 



If the system is in a non-equilibrium steady state then 
work should be done to keep it there. It is just the work 
needed for creating macroscopic currents inside the sys- 
tem. To illustrate this thesis, we can employ relation 
(0) for constant Ti, T 2 , divide it by dt and write it as 
F = W — II. The positive quantity II is the energy dis- 
sipated per unit of time. Using Eq. (|9j) we get 

F = (Ti - T 2 ) J d Xl dx 2 J 2 <9 2 lnP(xi|a; 2 ) - 



2 1 

J dx 1 dx 2 P(x 1 ,x 2 ) Y.( dlH + TA lnP)S 



(16) 



The last term in the right-hand side is nothing else but 
the energy dissipated per unit of time; it can be written 
alternatively as the sum of energy dissipation driven by 
the corresponding thermal baths: II — Y^i=i(Ti&iS — 
diQ)/dt (recall that diQ is the heat obtained from the 
thermal bath i, and &iS is the change of system's en- 
tropy induced by this thermal bath). The first term in 
the right-hand side of Eq. (|l^) should be associated with 
the performed work. In the stationary state the free en- 
ergy is constant, and the dissipated energy and the per- 
formed work are equal. 
Using Eqs. (Opl) we get 



Stat = 7 7^([^2] 2 )l " 7 2 r r {SF 2 d 2 A) , 
i 2 l 1 1 1 



IT K /L?ZT 12\ 2«( 2T 2 - Ti) 

n = 7jr(F2] )i-7 ^ 



(SF2d 2 A) , (17) 



where k — (Ti — T 2 )/Ti, and (...)o(i) means averaging by 
the distribution Tb(i)- We observe the following deceiv- 
ingly simple relation, valid to leading order in 7, 



n = T 2 s tot + o( 7 2 



(18) 



For a usual non-stationary system tending to equilibrium 
we have the following relation between entropy produc- 
tion and energy dissipation: II = TStot, where T is the 
temperature of the unique thermal bath. On the other 
hand, Eq. (^8|) reflects degradation of the energy in the 
stationary state. The distinguished role of T 2 is con- 
nected with conservation of detailed balance for small 
time-scales (see Eq.(||)). Indeed, both entropy produc- 
tion and energy dissipation are small on the characteristic 
times of the a; 1 -variable (|l7|). This equation also shows 
that when T 2 is close to zero, the energy dissipation (but 
not the entropy production) looses its leading term. At 
least in this limit the 7 2 -correction to Stot is negative. 

Let us apply the obtained general results to a sim- 
ple toy model. We consider a pair of weakly-interacting 
oscillators with coordinates x±, x 2 and Hamiltonian: 
H = oxi/2 + axl/2 + gx\x\, where a > 0, g > 0. Very 
similar models are applied to describe an oscillator with 
random frequency g or some electrical circuits jjj. For 
simplicity we shall discuss the model keeping only the 
first non- vanishing order in the small parameter g. The 
stationary distribution has the form (|l0|), with 



■(1 -ap 2 x 2 2 )(l -a(3ix\) 



g{Ti - T 2 



After some calculations we get from (O 
• 8g 2 (Ti - T 2 f 2 



(19) 



(20) 



n = % -f [Tl 3 T2) \ lT2 + 7 2 (T! - 2T 2 )) (21) 
1 1 a 6 

We see that in this model the 7 2 correction to Stot is 
always negative. For II it is only the case if Ti < 2T 2 . 
We have thus provided a concrete example of our general 
results. 

(iv) Let us now discuss the Onsager relations concern- 
ing heat transfer in the glassy state. These fundamental 
and experimentally testable relations were proposed by 
Onsager to describe transport in weakly non-equilibrium 
systems (the linear case) |ij [Q. Later they were gen- 
eralized to the non-linear regime. Following standard 
arguments W [|l0| the Onsager relation reads in our case 



3 



(22) 



where Qi, given by Eq. (|J), but see also Eq. (|l|), is the 
heat flux from the thermal bath i. In the stationary case 
one has Qi + Q 2 — 0. For our toy model we get from 
Eqs. (BUIT 



7 



g 2 ft -ft 
r ia 3 ft?/3 2 2 



(23) 



The linear case corresponds to Eq. ( p2| ) with ft ~ ft- 
Indeed, then the fluxes can be written in more familiar 
form: Qi = J2j LijAfij, where Aft = ft — ft is a small 
deviation of the inverse temperature ft from its equilib- 
rium value ft, and the Lij depend only on ft but not on 
ft. 2 separately. In that case the relation (^2|) takes the 
form: L12 = L 2 i (= 8g 2 Tg/r 2 a 3 in our toy model). Let 
us stress that this form of Onsager relations is applicable 
only for the linear case JjJ (!(]], in contrast to the more 
general relation (p2|). 

In fact the general validity of Eq. ( p2| ) in the linear 
regime is a fundamental theorem jlj supported by very 
general arguments. It means that any breaking of Eq. 
( p2| ) can be connected only with T\ ^ T 2 . The converse 
is not true: there are physically important cases when the 
Onsager relations hold for Ti ^ T 2 M. Thus checking 
these relations for our concrete class of non-equilibrium 
systems seems very important. From Eqs. (p"l]), (|T2|), 
(U) we obtain to leading order in 7 

dp x Q 2 - dfoQi = 

ft -ft 



experimentally. In glasses T 2 will correspond to effec- 
tively generated temperature 0- ||, and T\ is the tem- 
perature of the environment. By changing the cooling 
rate, also for these systems Eq. ( p2] ) constitutes a rela- 
tion between measurable quantities. A breaking of the 
Onsager relations can be investigated also in this more 
realistic setup. The details will be reported elsewhere. 

In conclusion, we have considered a stochastic model 
which contains all necessary ingredients of steady glassy 
behavior. In spite of the fact that such a system can 
be very far from equilibrium, it allows a thermodynamic 
description. Generalizing the usual equilibrium thermo- 
dynamics, we show that the glassy one can be obtained 
by minimization of the energy, keeping all entropies fixed. 
Universal relations, Eqs. (17, Il8|), are obtained between 



entropy production and energy dissipation. Energy dissi- 
pation (in contrast to entropy production) is closely con- 
nected to the fluctuations of slow degrees of freedom, and 
it looses its leading term when the corresponding temper- 
ature is close to zero. We discuss cases where the correc- 
tions arising from a finite but small ratio of characteristic 
times lead to a decrease of the entropy production and/or 
energy dissipation. We show that the nonlinear Onsager 
relation for heat transfer in the steady glassy state is al- 
ways broken, reflecting its strongly non-equilibrium char- 
acter. As the effect is of order unity, this breaking should 
be testable experimentally. It is also reminiscent of the 
breakdown of the Maxwell and Ehrenfest relations in the 
glassy state § §. 

A.E. A. is grateful to FOM (The Netherlands) for fi- 
nancial support. 
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{d Pl (T 2 ((5F 2 y)) +d P2 {T 2 ({5F 2 Y))} (24) 



In the linear regime with ft rj ft the Onsager relation 
is satisfied automatically. However, for the considered 
glassy system it is the exceptional case, and ( p2|) cannot 
be true in general. Indeed, for our toy model we get 



dp 1 Q 2 -d 02 Q 1 = 1 ^T 1 T 2 {T* 



T4 



(25) 



implying a violation of the Onsager relation for any 
Ti 7^ T 2 , because, due to (|2^), the right hand side of 
( p5| ) has the same order of magnitude as the individual 
terms in the left hand side. In this sense the violation 
is strong. One can give also general, model-independent 
arguments supporting breakdown of (|22|). Indeed, if it 
were to be valid, Eq. (^J) says that we should have 
{{6F 2 ) 2 } = ft /(ft - 02) for all U 2 , where / is some 
positive function. Such a form cannot hold for a triv- 
ial reason: taking the limit (3 2 — » one obtains zero in 
the right-hand side, while the left-hand side typically di- 
verges, or at least stays finite and non-zero including non- 
typical cases. Our prediction for the breakdown of the 
Onsager relations in the glassy state should be testable 
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